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We present a strong-coupling expansion of the Bose-Hubbard model which describes both the 
superfluid and the Mott phases of ultracold bosonic atoms in an optical lattice. By performing 
two successive Hubbard-Stratonovich transformations of the intersite hopping term, we derive an 
effective action which provides a suitable starting point to study the strong-coupling limit of the 
Bose-Hubbard model. This action can be analyzed by taking into account Gaussian fluctuations 
about the mean-field approximation as in the Bogoliubov theory of the weakly interacting Bose 
gas. In the Mott phase, we reproduce results of previous mean-field theories and also calculate the 
momentum distribution function. In the superfluid phase, we find a gapless spectrum and compare 
our results with the Bogoliubov theory. 
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I. INTRODUCTION 

Recent experiments on ultracold trapped atomic gases 
have opened a new window onto the phases of quan- 
tum matterfA^, A gas of bosonic atoms in an optical or 
magnetic trap has been reversibly tuned between super- 
fluid (SF) and insulating ground states by varying the 
strength of a periodic potential produced by standing 
optical waves. This transition has been explained on the 
basis of the Bose-Hubbard model with on-site repulsive 
interactions and hopping between nearest neighboring 
sites of the latticed As long as the atom-atom interac- 
tions are small compared to the hopping amplitude, the 
ground state remains superfluid. In the opposite limit 
of a strong lattice potential, the interaction energy dom- 
inates and the ground state is a Mott insulator (MI) when 
the density is commensurate, with an integer number of 
atoms localized at each lattice site. 

The Gross-Pitaevskii equation or the Bogoliubov 
theory 4 assume quantum fluctuations to be small and 
are unable to describe the SF-MI transition and the MI 
phase. The SF-MI transition is usually studied within a 
strong-coupling perturbation theory which assumes the 
kinetic energy to be small and treats exactly the on- 
site repulsion. In the simplest version, the kinetic en- 
ergy term is considered within mean-field theory. 3,5 ' 6,7 
The mean-field approximation is well known to give a 
reasonable estimate of the critical on-site repulsion at 
which the MI-SF transition occurs. Fluctuation correc- 
tions to the mean-field approach have also been consid- 
ered within a systematic strong-coupling expansion^ All 
these approaches have given a reasonable description of 
the MI phase and in particular of the excitation spec- 



trum. However, they have not provided a description of 
the SF phase. 



In this work, we develop a strong-coupling expansion 
of the Bose-Hubbard model which allows us to extend 
the treatment of Refs. and describe both the 

MI and SF phases. Our approach is similar to strong- 
coupling expansions introduced for the (fermionic) Hub- 
bard modelA^i In Sec.[H] we derive an effective action for 
the Bose-Hubbard model in the strong-coupling limit by 
performing two successive Hubbard-Stratonovich trans- 
formations of the intersite hopping term. This effective 
action involves the exact one- and two-particle Green's 
functions in the local limit (i.e. in the absence of inter- 
site hopping). We then use the standard Bogoliubov ap- 
proximation: we perform a saddle-point (or mean-field) 
approximation and expand the action to quadratic order 
in the fluctuations (Sec. lIIIJl . In the MI phase, we recover 
the previous mean-field result i*^ We find a gapped ex- 
citation spectrum which becomes gapless at the MI-SF 
transition. We also calculate the momentum distribution 
function and study the critical behavior at the transition. 
In the SF phase, we obtain a gapless spectrum (in agree- 
ment with Goldstone theorem) and compute the Bogoli- 
ubov sound mode velocity. We compare our results with 
the Bogoliubov theory. 
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II. EFFECTIVE ACTION IN THE 
STRONG-COUPLING LIMIT 

The Bose-Hubbard model is defined by the Hamilto- 
nian 



H 



t E + h.c.) -/x^ri r + ^^n r (n r - 1), 

(r,r') r r 

w 

where ^„^J are bosonic operators and n r = ^J^ r - The 
discrete variable r labels the different sites (i.e. minima) 
of the optical lattice, t is the hopping amplitude between 
nearest sites (r, r') and U the on-site repulsion. The op- 
tical lattice is assumed to be bipartite with coordination 
number z. The density, i.e. the average number n of 
bosons per site, is fixed by the chemical potential fi. 

We write the partition function Z as a functional in- 
tegral over a complex field if> with the action S[ip*,ip] = 

dr{J2 T V'^rV'r + H[ip*,ip}} [t is an imaginary time 
and P — l/T the inverse temperature]. Introducing an 
auxialiary field <j) to decouple the intersite hopping term 
by means of a Hubbard-Stratonovich transformation)*^ 
we obtain 



Z / V 



-(4,\t- 1 <f')+W[r,4>] 



(2) 



where we use the shorthand notation (<f>\ip) = J2 a 4>*a'4 > a. = 
Jq dr a J2 r 4>* ( v a)ip( v a)- t 1 denotes the inverse of the 
intersite hopping matrix defined by t rr * = t if r, r' are 
nearest neighbors and t rr > = otherwise. So and Zq 
are the action and partition function in the local limit 
(t = 0). (•••)o means that the average is taken with 
Sq[4>*, ip]. In the last line of (J5J, we have introduced the 
generating function = ln(exp J^ai^a^a + c.c.))o 

of connected local Green's functions*^ 
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where {a,i,bi} = {ai 
we obtain 
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W[<f>*,<f>] = E 

R=l 

where means that all the fields share the same value 
of the site index. If we truncate W[(j>* , <fi] to quartic order 



in the fields, we obtain the action 



a . b 
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where G = G 1 . Eq. Q was used as a starting point by 
van Oosten et. al. to study the instability of the MI 
with respect to superfluidity^ Their results are summa- 
rized in Appendix El and lead to the usual mean-field 
phase diagram shown in Fig. It is tempting to go be- 
yond the mean-field approximation by considering Gaus- 
sian fluctuations of the <fr field about its mean-field value. 
The Green's function obtained in this way is however not 
physical since it leads in the SF phase to a spectral func- 
tion which is not normalized to unity^ Physical quan- 
tities like the excitation spectrum, the velocity of the 
Bogoliubov sound mode or the momentum distribution 
in the SF phase are therefore out of reach within this 
approach. 




FIG. 1: Phase diagram of the Bose-Hubbard model showing 
the superfluid phase (SF) and the Mott insulating (MI) phases 
at commensurate rilling n. The dashed lines corresponds to a 
fixed density n = 0.2, n = 1 and n = 2. For a commensurate 
density n, the MI-SF transition occurs for U/(zt) — 2n + 
1 + 2(n 2 + n) 1/2 (for n = 1, this yields U/{zt) ~ 5.83, i.e. 
U/t ~ 23.31 for a two-dimensional atomic gas in a square 
optical lattice). 

These difficulties can be circumvented if one performs a 
second Hubbard-Stratonovich decoupling of the hopping 
term: 



• Or, &i • • • £>h}. Inverting Eq. 10, g 



Z / V[if>*, <f>]e&W)-l(*W^-]+m4>*,4>]. (6) 



In Appendix El we show that the auxiliary field of this 
transformation has the same correlation functions as the 
original boson field (hence the same notation for both 
fields). The effective action S[ip*,ip] is obtained by inte- 
grating out the (j) field in Eq. © . This procedure was car- 
ried out in detail in Ref. in the context of the fermionic 
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Hubbard model. Similarly, we obtain^ 

a , b 

+ \ ^i^MM^ai^^^b,, (7) 

ai ,02,61,62 

where T u (ti, T2\ T3, T4) is the (exact) two-particle vertex 
in the local limit. In Eq. (JJ, we have neglected R- 
particle vertices (R > 3) whose amplitudes are given by 
the (exact) local i?-particle vertices r R iAfi T 11 is local in 
space but has a complicated time dependence (see Ap- 
pendix [Bj • In the following, we approximate T n by its 
static value (obtained by passing to frequency space and 
putting all Matsubara frequencies to zero) . This approx- 
imation is justified for energies much below U where the 
frequency dependence of the local two-particle vertex is 
weak. At higher energies, its validity is more difficult to 
assess. Introducing 



(8) 



we finally obtain 



S = - 



f dTdT'Y,rAr)[G- 1 {v,r;v',T') 
a f p 

+U,r>5(T - r')] W (r') + | / dr V V; *P*r*PrA- 

z Jo r 

(9) 

The action is the starting point of our analysis. 

It is analog to the original action J Q dT{^2 T ip*d T ip T + 
H [ip*, ip]} with two noteworthy differences: the "free" 
propagator involves the exact local propagator G, and 
the amplitude of the boson-boson interaction is given by 
the exact local two-particle vertex (approximated here 
by its static limit) . The action yields the exact parti- 
tion function Z = Zq J T>[ip* , ip]e~ s and the exact Green 
function — {ip r {T)ip*,{T')) both in the local (t — 0) and 



non-interacting (U — 0) limitsi2*AS. By means of two 
successive Hubbard-Stratonovich transformations of the 
intersite hopping term, we have thus performed a partial 
resummation of interaction processes and obtained an ef- 
fective action which provides a suitable starting point in 
the strong-coupling limit. 



III. MEAN-FIELD AND GAUSSIAN 
APPROXIMATIONS 

In order to study the Mott and superfluid phases 
from the strong-coupling effective action 10, we use the 
standard Bogoliubov approximation: we first perform a 
saddle-point (or mean-field) approximation and then ex- 
pand the action © to quadratic order in the fluctuations. 
The saddle-point action is given by 



^ = -(G- 1 +^ 2 + f< 



(10) 



where G = G(iuj = 0), D = zt, and N is the total number 
of lattice sites. The saddle-point value ipo (assumed here, 
with no loss of generality, to be real) is obtained from 

dS/dip = 0: 



G~ 



D 



if G~ 



D > 0, 
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(11) 



otherwise. 



The MI-SF therefore occurs when G _1 + D = 0, in 
agreement with the results of Appendix El which leads 
to the phase diagram shown in Fig. Using (ip T ) = 
5lnZ(J*,J)/5J*\j*=j=o, where Z[J*,J] is given by 
Eq. IjAljl of appendix El we obtain (f>o — Dipo where 
0o is the mean value of the auxiliary field. Near the 
MI-SF transition, where G -1 + D w 0, we then find 
0o — 2(Z?~ 1 + G)/G IIc in agreement with the result of 
Appendix [0 

To quadratic order in the fluctuations ipr — — ipo, 
we obtain the action 



54y(f(k,« W )i(-k,- W ))f'f^ +£k+2 * 2 ^ . u ^^nff'^ A (12) 

2/ -^rv, hYK ' "\91pQ -G 1 (-iw) + e_ k + 2gVo J \ ip* (-k, -iu) J K ' 

k,o; 

1 

where ^(k, iuj) is the Fourier transformed field of ip r (r) A. Mott phase and the MI-SF transition 

and oj a bosonic Matsubara frequency, ek, the Fourier 
transform of — t TtT i, is the boson dispersion in the absence 
of the one-site repulsion. 

In the Mott phase, where ipo — 0, the Green's func- 
tion Q(k, iuS) — — (ip(k, iui)ip*(k, iui)) can be directly read 
off from Eq. (Q2J: g- l (k,iuS) = G-\iu;) - e k . Using 
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Eq. JB2ll . one obtains 



1 - z k 

ilo — EC 



2k 



itu - Ei 



(13) 



The two excitation energies E k and the spectral weight 
Zk are defined by 



k p 2 2 



+ 4e k f/x + C/ 2 



1/2 



E£ + Sfi + Ux 



(14) 



where a; = no + 1/2 and (5n = (X—U (no — 1/2) . no e% (/x) 
is the (integer) number of bosons in the local limit for a 
chemical potential « (see Annendix lBjl . 

The excitation energies E£ , and the correspond- 
ing spectral weight Zk and 1 — Zk, are shown in Figs. [2} 
El in the MI n = 1 of a two-dimensional atomic gas in 
a square optical lattice. The spectrum exhibits a gap 
E+ =0 - E k=0 = (D 2 - ADUx + U 2 ) 1 ' 2 which decreases 
as U decreases. The MI becomes unstable against su- 
perfluidity when E k=Q = or E^_ = 0, which agrees 
with Eq. ilC3j) of Appendix [0 an d leads to the phase 
diagram shown in Fig. ^ The gap E ] 



k=0 



k=0 



{D 2 -ADUx+U 2 ) 1/2 closes at the transition if both E+ =Q 
and E^_ Q vanish, which occurs at the tip of the Mott 
lob. The MI-SF transition then takes place at fixed den- 
sity, which is the situation of physical interest. Figs. EJ 
are obtained with a chemical potential = —D/2, 
which ensures that the MI-SF transition takes place at 
fixed density n = 1 (see Appendix E|. The decreasing 
of the Mott gap is accompanied by an increase of spec- 
tral weight at k = 0, which diverges at the transition. 
Figs. EE also show the results of the Bogoliubov theory 
(as applied to the original Hamiltonian QJ). The Bo- 
goliubov theory always predicts the ground-state to be 
superfluidifi Away from k = 0, it provides a good ap- 
proximation of the negative energy branch E^ but gives 
a poor description of E£ . 

If we expand the equation E k=0 — to order 0{t 2 /U), 
we obtain 

D 2 

u - C/n + L>(n + 1) + — (n 2 +n ) = 0, 
D 2 

H-U{n -1)-Dn -— (n 2 +n ) - 0, (15) 

which differs from the energy calculation of Ref. |g by 
terms of order 0(t 2 /U). This discrepancy results from 
the neglect of the one-loop correction due to T 11 in the 
calculation of the Green's function [Eq. lTL3"l) ]. which also 
gives a contribution of order 0(t 2 /U). However, even 
without this term the phase diagram looks qualitatively 
similar to the Freericks and Monien phase diagram. 

From the Green's function lllot . we can also ob- 
tain the momentum distribution nk = (ip^/ipk) = 

— J° dwyl(k, uj) = 1 — Zk- nk measures the spectral 
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FIG. 2: Top: Excitation energies E£ (solid line) and 
(dashed line) in the MI n = 1 for U = 30t. Bottom: Spectral 
weight Zk (solid line) and 1— Zk (dashed line). The dotted lines 
show the result obtained from the Bogoliubov theory (which 
predicts the phase to be superfluid). [F = (0,0), M = (tv, it) 
and X = (7r,0).] Results shown in Figs. 12151 are obtained for 
a two-dimensional atomic gas in a square optical lattice. 
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FIG. 3: Same as Fig. [2 but for U = 25i. 



weight of the negative energy E^ of the spectrum. Deep 
in the Mott phase, the momentum distribution is roughly 
flat. Closer to the MI-SF transition, a peak develops 
around k = 0. This peak diverges at the transition 

(Fig. a. 

The critical theory of the SF-MI transition can be ob- 
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Momentum distribution 




FIG. 4: Momentum distribution n k = (V'k^/'k) in the MI n — 1 
for U = 30t (top) and U = 25 (bottom). 



tained from the action ^ by expanding the inverse prop- 
agator G~ 1 (io;)— ek to quadratic order in k and to. Noting 
that dG (ito) /d(ito)\i u= o = dG~ x jd[i (and similarly for 
the second-order derivative), we obtain 



S 



where 



dr \ dr 



+K 3 \V^ T \ 2 + -\^ r \ 4 



ro oc G 



T-olV^r I + K^rdrA + K a \d T i^\ 2 
u , 



(16) 



D, 



Ki oc — . 

OjJb 



(17) 



At all points on the MI-SF transition line except at the 
Mott lob tip, ro vanishes but K\ remains finite. The 
critical theory has then a dynamical exponent z = 2. At 
the tip of the Mott lob where both r$ and K i vanish, the 
dynamical exponent z = 1. A similar analysis, based on 
the effective action 5 [</>*, 0], can be found in Ref. Q. 



B. Superfluid phase 

In the SF phase (tpo 7^ 0), the Green's function of 
the ip field is obtained by inverting the 2x2 matrix 
propagator in Eq. Ill2t . For the diagonal component 



£/(k, ito) = 
Q{k,ito) 
where 

e£ 2 = 

A k = 
S k = 

S k = 
G k = 



■(^(k, iu))ip*(]t, ito)), we obtain 

(ito + 8/i + Ux)(ico — z£)(iio — z k ) 
= (^+E+ 2 )(u 2 + E^ 2 ) : 



(18) 



4C k ) 1/2 , 



^±i(i 2 -4B k )^, 
2^-2(G- 1 +D)-e k , 
-(2G" 1 + 2D + e k )(fy + £/x) 
(G^+Z}) 2 , 



r/ 2 



2S k 

£? k - (G _1 + D) 2 (S^l + Ux) 2 



(19) 



From ltl8)l . we deduce the spectral function A(k, to) 



A(k,w) = 

(£+ + <5// + E/:e)(£;+ 



W^k) 
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2£+(£+ 2 -£ ] 
^-l^)(£+ + z+)(£+ + z k ) 
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+ 8n+Ux){E£ 


"4)(^k - 




2E~(E+ 2 




(£k 


-8(i-Ux)(E k 


+ 4)(i? k - + z k ) 




2E U (E+ 2 


-£ k 2 ) 



<y(w + e+ 



5(uj 
8 (to 



(20) 



The Green's function l|18|) has the desired physi- 
cal properties. The spectral function is normal- 
ized, X^o^ w A(k, to) = 1, and has the correct sign: 
sgn[A(k, to)] = sgn(w)il2i There are four excitation 
branches ±-E k , two of which (i-EjJT) being gapless for 
k — > (Fig. EJ. However, for a given value of k, only 
two branches carry a significant spectral weight. Away 
from k = 0, the spectral weight is almost completely ex- 
hausted by E£ and — E^. In the vicinity of k = 0, the 
two gapless branches ±-E k exhaust the spectral weight. 
By expanding E^ in the vicinity of k — » 0, we find a 
linear spectrum 



E: 



clkl 



(21) 



where 



2t{G- 1 + £>) 



1/2 



: 2 + 2 7 (G- 1 +£) 

S^i 2 + 2SfiUx + U 2 /i 
(Sn + Ux) 2 ' 
U 2 (x 2 -l/4) 
(8(i + Ux) 3 ' 



(22) 
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PIG. 5: Excitation energies ±E^ and spectral weight in the 
SF phase n = 1 and U = 20. 



Our strong-coupling approach therefore reproduces the 
Bogoliubov (Goldstone) mode of the SF phase. 

As discussed in Sec. IIII Al our strong-coupling theory 
is not an expansion order by order in t/U, For this rea- 
son, the computation of the chemical potential from the 
single-particle Green's function, i.e. n = Tr(Q), is not 
reliable. We have therefore used the chemical potential 
obtained within the mean-field approximation discussed 
in Appendix El 

Fig. also shows the results of the Bogoliubov theory 
(as applied to the Hamiltonian Q) for the same chemical 
potential p. The Bogoliubov theory provides a good ap- 
proximation to E^ and therefore to the low-energy part 
of the excitation spectrum. This implies that the veloc- 
ity of the gapless mode [Eq. ll22l) ] can be approximated 
by the Bogoliubov result c = [2t(p + D)] 1 / 2 . Away from 
k = 0, the Bogoliubov approach gives a rather poor de- 
scription of . 

The Green's function Q(k, iuj) yields the momentum 
distribution 

rik = ($kVk) 

,o 

= Nifa&kja - / doA(k,w), (23) 

J — oo 

Apart from the condensate contribution Ntp^Skfi, the 
momentum distribution function is directly given by the 
spectral weight of the negatives energies -E+ and -££ 
(Fig.EJ. 

Fig. shows the integrated spectral function p(uj) = 
I (27r)'a v) for a commensurate density n = 1. Deep 
in the Mott phase, p(uj) is essentially given by the non- 
interaction density of states of free bosons on the square 



lattice centered around —fj, and U — p and with relative 
spectral weigths —no and no + 1. The two peaks near lo = 
—p and u> = U—fj, are due to the Van Hove singularities in 
the density of states of free bosons. When decreasing the 
value of U/t, the Mott gap decreases and p(u>) strongly 
increases at the gap edges. At the critical value U/t ~ 
23.31, the gap closes and p(ui) diverges at u> — 0. This 
divergence persists in the superfluid phase. 



IV. CONCLUSION 

By performing two successive Hubbard-Stratonovich 
transformations of the intersite hopping term, we have 
shown how to derive an effective action which provides a 
suitable starting point to study the strong-coupling limit 
of the Bose-Hubbard model. This action can then be 
analyzed by taking into account Gaussian fluctuations 
about the mean-field approximation as in the Bogoliubov 
theory of the weakly interacting Bose gas. The main im- 
provement over previous related approaches£»£*£»& is the 
possibility to describe both the Mott and SF phases. 
Both in the Mott and SF phases, we compute the excita- 
tion spectrum and the momentum distribution. Our ap- 
proach clearly shows how the excitation spectrum, which 
is gapped in the MI phase, becomes gapless at the MFSF 
transition. 

The strong-coupling expansion presented in this paper 
should in principle also applies to more complicated situ- 
ations where for instance several atom species are present 
in the optical lattice. 

Note added: after completing this paper, we became 
aware of two related works. Konabe et al. 14 have stud- 
ied the single-particle excitation spectrum in the Mott 
phase and obtained results similar to ours. The method 
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used by these authors bears some similarities with the 
strong-coupling expansion discussed in the present pa- 
per. Within a slave-boson representation of the Bose- 
Hubbard model, Dickerscheid et. al*&> have discussed 
both the Mott and SF phases. Their results agree with 
ours (whenever the comparison is possible). 



where J*, J r are external sources. After the Hubbard- 
Stratonovich decoupling of the intersite hopping term [see 
Eq. 0] and the shift <j>* -> <p* - J* ,<j> -> <j> - J of the 
auxiliary field, we obtain 



APPENDIX A: HUBBARD-STRATONOVICH 
TRANSFORMATIONS 

The Green's functions of the boson field ip can be ob- 
tained from the generating function^ 

Z[J\J] = J V[4>\^]e Wt ^- So[r ^ ]+[{JW+c - c \ (Al) 



J 



Z[J\J] = J D[^*,^,^*,0] e -(*- J l*~ 1 C*- , ))+K*W+ c - c -]-' sr «»W.*] 

= Z [ T>[<p*,4>]e-^- J ^ 1 ^- J ^ +w ^^. (A2) 



A second Hubbard-Stratonovich decoupling of the hopping term (with an auxiliary field ip') leads to 



Z[J*,J] = ZojD^V^',*'*'!'*''-^^-^^^^^' 

= Zq [ P[vyV>*,0]e W ' l * V '' ) - [( ' A ' l0)+c - c - 1+[( ' / '' |J)+c - c - 1+H ' [0 *' 01 . (A3) 



From l|A3Jl we deduce that Z[J* , J] is also the generating 
function of the Green's functions of the ip' field, ip' can 
therefore be identified with the original boson field ip. 



APPENDIX B: CALCULATION OF THE LOCAL 
GREEN'S FUNCTIONS G AND G u 

In the absence of intersite hopping (t = 0), the states 
\p) = (p > integer) are eigenstates 

with eigenvalues e p = —fip+ (U/2)p(p — 1). [We consider 
a single site and therefore drop the site index.] |0) is the 
vacuum of particles. This yields the partition function 



Zn 



p=0 1 



-0ep 



. In the ground-state, for a given value 



of the chemical potential fi, there are no bosons per site, 
where rig is obtained from e„ = min p e p . The latter 
condition leads to rto — 1 < Ia/U < hq if /i > —U, and 
no = if /i < — U. Note that n is integer (except when 



fi/U = p is integer; the states \p) and \p + 1) are then 
degenerate), even when the boson density n is not. 

The single-particle Green's function G(r) = 
— (Trip^)^ (0)) is easily calculated using the clo- 



relation J2^Lo\p)(p\ = 1- F° r T > 0, 



finds 



G(r)^~—Y / (P+^- ( ^ T)ep - TEp+1 , (Bl) 



u p=0 

and, in frequency space, 
-n 



n Q + 1 



iuj + e no _i — e„ iu + e no — e na +i 



(B2) 



where to is a bosonic Matsubara frequency. 

The two-particle Green's function can be calculated in 
the same way. One finds 
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G II (r 1 ,T 2 ;T3,r 4 = 0) = (T r ^(Ti)^(r 2 )^(0)^t( r3 )) 



-0e p 



(p +l)(p + 2)e Tl ^- e p +l)+T2(€p + 1 - €p + 2)+T3{e * >+2 - e *'+ l) 9(T 1 - t 2 )0(t 2 - t 3 ) 



p=0 



+{p + l)(p + 2)e Tl(£ p+ 1 -^+ 2)+r2(e ''- c ''+ l)+r3(e p+ 2 - e ''+ l) 6'(T 2 - 71)0(71 - r 3 ) 

+ (p + l)2 e rx(e p -e p+1 )+r 2 (e p -e p+1 )+T3(^ +1 -e p )^ (n _ _ ^ + ^ _ _ n) ] 

+p(p + i)e T ^-^^+ T2 ^^ +T ^~^-^9{T^ - 7-i)fl(Tx - r 2 ) 

+ i) e ^(^-^+i)+^( ep -i- ep )+r 3 ( ep - ep - 1 ) 6 ,( r3 _ T2 )<9( r2 _ n )l . (B3) 

I 



After a somewhat tedious calculation, we obtain for the 
Fourier transform of the connected part in the static 
limit: 







4(n + l)(n + 2) 


(2/i- 


-(2no + l)C0(^"o- 


-m) 2 




4n (n - 1) 




(/*- 


E/(n - l)) 2 (*7(2n 


-3) 




4n (n + 1) 




(m- 


Un ){-fi + U(n - 


I)) 2 




4n (n + 1) 




0- 


[/n ) 2 (-/i + £/(n - 


-1)) 




4n 2 




(-/* 


+ C/(n -l)) 3 




4(n 


+ 1) 2 




(a*- 


C/no) 3 ' 





(B4) 



The static limit of the two-particle vertex T n is equal to 
-G Uc /G 4 . 



(see Appendix El . The ground-state energy per site E 
— lim^oo -J-g In Z is then given by [see Eq. 0] 



dT 1 dT 2 dT 3 G u (T 1 ,T 2 ;T 3 ,0) - 2f3[G(iu = 0)] 2 where a 



E = a + a 2 (j>o + 04^0, (C2) 
- lim^^oo In Zq is the ground-state en- 
ergy in the local limit, a 2 = D^ 1 + G, and 04 = — jG Uc . 
The mean-field value </>o is obtained by minimizing E. <j)Q 
vanishes in the Mott phase (a 2 > 0) and takes a finite 
value in the SF phase (a 2 < 0). The MI-SF transition is 
then given by a 2 = 0, which leads to 



D 1 

"2 ±2 



D +U — 4DUx 



1/2 



(C3) 



where no is the integer number of bosons in the local 
limit for a chemical potential /t (see Appendix [5J| ■ x 
and Sfi are defined in Sec. IIIII For each value of no, 
Eq. l|C3Jl defines a Mott lob in the U — /1 phase diagram 
(Fig. nj, whose tip corresponds to <5/t+ = <5/i_ = —zt/2 
and U/(zt) = 2n + 1 + 2(n§ + no) 1 / 2 . At the lob tip, 
da 2 /d/i = 0. 

In the SF phase, the order parameter <j)Q is given by 
0q = —02/(204), and the ground-state energy takes the 
value 



APPENDIX C: AUXILIARY-FIELD MEAN-FIELD 
APPROACH 

In this appendix, we review the mean-field results 
obtained from the action S[4>*,<j)} [Eq. iJSJ]- 6 Within a 
saddle-point approximation, where the field <fio is taken 
real and assumed to be time and space independent, we 
action becomes 



S 

N/3 



{D- l +G)4>l 



1 



:G 



(CI) 



where D — zt. G and G IIc are the single-particle and 
two-particle local Green's functions in the static limit 



E 



a 



4a 4 



From (£31), we deduce the mean boson density 



dE 



1 d 



= »o 



a 2 da 2 



(C4) 



4 dfi ^ 04 ) ° 2a4 <9/t ' 



where the last equality holds near the MI-SF transition 
(a 2 « 0). We have used no = —dao/dfj,. We conclude 
that, at the MI-SF transition, the boson density remains 
pinned at the integer value no if da 2 /dfi = 0, which cor- 
responds to the tip of the Mott lob in the /t — U phase 
diagram (Fig.QJ. 
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